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Abstract: Recently, the Fokker-Planck dynamics of particles in periodic 
potentials ±V, have been investigated by using the matrix continued fraction 
method. It was found that the two periodic potentials, one being bistable and the 
other metastable give the same diffusion coefficient in the overdamped limit. We 
show that this result naturally follows from the fact that the considered potentials 
in the corresponding Schrodinger equation form supersymmetric partners. We 
show that these differing potentials ±V also exhibit symmetry in current and 
diffusion coefficients: J+{F) = — J_(— F) and D+(F) = D_(—F) in the presence 
of a constant applied force F. Moreover, we show numerically that the trans- 
port properties in these potentials are related even in the presence of oscillating drive. 
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I. INTRODUCTION 



Periodic potentials are common occurrences in diverse kinds of extended systems 
such as crystals in solids and microtubules in living systems. The nature of the 
potential profile depends on the medium and the particle moving through it. A Ag + 
ion, for example, of a Agl crystal in the superionic phase [1] moving in a certain 
direction may be gainfully modeled to encounter a periodic potential profile V(x) = 
A cos kx + B cos 2kx with two wavevector components instead of a simple cosinusoidal 
potential. Such potential profile modeling could be based on the consideration of 
crystal structure and the underlying energy landscape in the direction of motion of 
the ion. In some other direction, however, the periodic potential profile could 
be entirely different. The potential profile —V(x) is also periodic with the same 
periodicity as V(x) yet they have important differences. The potential V(x) is bistable 
in nature, whereas —V(x) has a metastable well within a period, Fig.l. 

Recent studies show that an overdamped Brownian particle exhibits identical dif- 
fusive behavior in the potentials V(x) and — V(x) in the entire range of parameter ^ 
[2]. However, the diffusion constant of underdamped Brownian particles are quite dif- 
ferent for the two related potentials. In this work we reason and explain that this result 
follows from the fact that the potentials of the Schrodinger equation with (Euclidian 
time) corresponding to Fokker-Planck equation of motion [3] with potentials ±V(x) 
form supersymmetric pairs and hence have identical eigenvalue spectra[4,5,6]. The 
identical eigenvalue spectra along with one to one correspondence between their eigen- 
functions ensure that both the systems (i.e., with ±V(x)) exhibit same long time 
average asymptotic behavior. This also implies that when a constant force F is applied 
to the potentials ±V(x) the average stationary current will satisfy J+(F) = — J-(— F) 
and the diffusion constant D + (F) = D_(—F). This result can also be seen from the 
analytic solution of the Fokker-Planck equation. In addition, we show that when the 
system is driven periodically, both the potentials show similar diffusive behavior. 



II. THE MODEL 



Consider the motion of a particle moving in a tilted periodic potential V(x) and 
described by the overdamped Langevin equation [3] 

where the Gaussian random force satisfies (£(£)) = and (£(£)£(£')) = 2^T5(t — t'). 
Here (...) denotes thermal average at temperature T(written for /cpT, where ks is the 
Boltzmann constant). 7 is the frictional drag coefficient. The coresponding Fokker- 
Planck equation (in dimensionless units) is written as [3] 

^^- = L FP P(x,t), (2) 

where the Fokker-Planck operator ~L F p is given by 

L , P = iL (1 / w+r iL). ,3) 

The prime over V(x) denotes derivative with respect to x. After constructing the 
Hermitian operator L = e* //2 Lppe~' I>//2 , with the effective potential <3>(x) = the 
time indepedent Schrodinger equation corresponding to the Fokker-Planck equation 
(2) can be written as [3] 

L^s) = [T— - V 8+ (x)\1> n {x) (4) 

Lip n (x) = -\ n ip n (x) (5) 

The Schrodinger eigenfunction ip n (x) = e~^ x ^ 2 (j) n (x), and the eigenvalues A n (> 0) are 
given by 

L FP (f) n ( X ) = -X n (j) n (x) (6) 

Here T substitutes for ^ in quantum mechanical Schrodinger equation. The 
Schrodinger potential Vs+(x) is obtained for the potential V(x) as 



For the potential — V(x) the corresponding Schrodinger potential is 

V" v' 2 

Vs-( X ) = T + W (8) 

The potentials Vs±(x) = W 2 (x)^fy/TW' (x) constitute the supersymmetric partner 
potentials and W(x) is referred to as the superpotential in supersymmetric quantum 
mechanics literature. The superpotential W(x) can also be expressed in terms of the 
ground state wavefunction of the Hermitian operator L [4]. For periodic potentials 
±V(x) in Fokker-Planck dynamics, it is known from supersymmetric quantum me- 
chanics that the eigenvalue spectra (with energy bands and band gaps) of Vs±(x) of 
the corresponding Schrodinger equation in quantum mechanics are strictly isospectral 
and the eigenfunctions are directly related to each other being in one to one corre- 
spondence [4]. The immediate implication of this observation is that the potentials 
±V(x) will exhibit identical transport behavior in the long time asymptotic limit. It 
is expected that the average or expectation values of all the dynamical variables for 
±V(x) will show similar behavior in the time asymptotic regime [4]. In fact for a given 
potential V(x) in the presence of applied field F the analytic expressions for current, J 
and diffussion coefficient, D, are given by[7,8] 

_ 1 - exp(- 27 rF/T) 

J * n rx +2w dxj / \ ' \ V ) 
Jxo 2tt±W 



and 



where, 



j* o+2w dxi±(x)i + (x)iJx) , s 

D=D ° ur* ^ ■ (10) 

1 V(x)-Fx f x -V(y) + Fy 

I+(x) = — e t / dye t ^ (11) 

L>Q Jx-2-K 
1 -V( x ) + Fx fX + 2-K V(y)-Fy 

I-(x) = — e T / dye T . (12) 
-Do Jx 

Here 2ir is the length of the period and D is the bare diffusion constant (= fcT/7). By 
tedious algebraic manipulation (for example see the appendix of ref [7]) one can show 
that, for potentials ±V(x), current J and diffusion coefficient D, respectively, satisfy 



the relations J + (F) = -J_(-F) and D + (F) = D_(-F). The subscripts ± refer to 
the potentials ±V(x). The result obtained in references [2] follows as a special case 
whenF = 0, namely D + (F = 0) = D_(F = 0). 

III. RESULTS AND DISCUSSION 

We have considered the potenitial V(x) = Acoskx + B cos2kx in the Langevin 
equation in equation(l) in its dimensionless form. We, in the following, as has been 
done from equation(2) onwards, measure energy in units of A, wavenumber in units 
of k and the friction coefficient in units of 7 and all other variables are given in terms 
of combinations of A, k and 7. Thus in these dimensionless units the potential reads 
as 

V(x) = cos a; + B cos 2a;, (13) 

etc. Thus a particular choice of B will fix the relative values of the two potential 
barriers . In a series of papers [2] , the ratio of the diffusion constant to the bare diffusion 
constant (-Do) has been given in the entire range (0,1) of the ratios of the potential 
barriers of V(x) in the overdamped as well as in the underdamped regime. In these 
works they calculate the diffusion coefficients by solving the Fokker-Planck equation (2) 
and (3) by the method of matrix continued fraction[3]. As mentioned in the previous 
sections their results in the high friction limit, can be understood from supersymmetric 
considerations as well as obtained from analytic expressions. We calculate in this work 
the diffusion constant D and current J in the presence of a constant force F. We 
calculate these quantities also for an asymmetric potential profile [fig. 2] 

V(x) = cosx + cos 2a; + 0.25 sin 3a;. (14) 

In the following the nature of currents J and the diffusion constant D calculated using 
expressions (9) and (10) respectively are presented. 

Fig. 3 shows the variation of current as a function of the constant applied force 
F at temperature T = 0.1 for the symmetric potential V(x) = cosx + cos 2a;. As 



expected, for low values of |F| << F c (where the largest barrier to motion just 
disappears), the current effectively remains zero as the largest barrier height is much 
larger than T = 0.1 (at F — the largest barrier height is 3.0). It may be noted that 
at temperature T = 0, there will be no current up to the critical field .F c (~ 2.74) as 
the particle will be in a locked state. Beyond F c the particle being in a running state 
leads to a finite current [3]. At finite temperature current begins to flow even before F c 
due to thermal activation. J approaches its linear regime (mobility /i — 1) soon after 
around F equals 3. From the figure it is clear that J±(F) = —J±(—F) as also J+(F) = 
J-(F), where the subscripts (±) on J corresponds to the potentials ±V(x). Our results 
also corroborate the earlier obtained result that the diffusion constant acquires giant 
enhancement [7,8] at the critical field F = F c [inset of fig. 3]. The giant enhancement 
in D near F c is expected as a fall out of the instability between the locked state 
(F < F c ) and running state[7,8]. We also observe that both the potentials ±V(x) 
exhibit identical diffussive behavior D + (F) = D-(F) for this symmetrical potential. In 
fig.4, we plot J± and D± for the periodic potential ±V(x), with V(x) = cosrc+cos 2x+ 
0.25 sin 3x, for the same temperature T = 0.1. Note that the crictical field now have 
different values (F cl fa —3.19 and F c2 ~ 2.12) for different applied field directions. With 
the asymmetric potential we see that the relation J + (+F) = — J_(— F) is satisfied as 
mentioned in the introduction. J(F) shows very similar behavior to what is shown 
in fig.3. The figure in the inset shows D + (F) = D-(—F). For both the potentials the 
diffusion coefficient shows the resonant enhancement close to the corresponding critical 
fields F c i and F&. And the diffusion coefficient does show the symmetry with respect to 
±F, for either of the potentials ±V(x): D(F) ^ D(-F), however, D + (F) = D_(-F). 

So far we have considered the case for particle subjected to constant force. We 
now study particle subjected to ac forcing where particle probability density does not 
evolve to a steady state(rather it evolves to a time periodic state). It may be noted that 
supersymmetric arguments do not apply when time dependent forcings are included. 
However, even in this case we show numerically that diffusion constant shows the same 



behaviour for potentials ±V(x). In Figs. 5 and 6 we plot the diffusion coefficient, as a 
result of an ac drive F = A cos(wt) , as a function of A at temperature T = 0.25. This 
is shown by solving the corresponding Langevin equation using Huen's method[9]. The 
diffusion constant D is obtained in the long time asymptotic limit as 

D=liml[(x 2 (t)>-(x(0> 2 ] (15) 

t^oo 2,1 

Where (...) denotes the ensemble averaging over a large number of trajectories. We 
discard the initial transients (t = 500r) and allowed the system to evolve for t = 
25000r. The time step r is taken equal to 0.01 and 3500 ensembles are used to evaluate 
the averages (15). 

In fig. 5, we give the plot for the symmetric potential V(x) = cosx + cos 2x for 
two frequencies ou — 1.0 and u = lO.O(inset). It is very clear from the figure that 
D is same for both the potentials ±V(x) within our numerical accuracy. Also, for 
uo = 1.0, D shows a resonant behaviour at 3.0 < A < 4.0. It should again be 
noted that the largest barrier height shown by the potential V(x) equals 3.0. That 
is, peaking behavior is exhibited close to the extreme forcing equal to the critical 
static field F c . This is where the particle has maximum likelyhood of diffusing forward 
as well as in the backward direction thereby enhancing the diffusion coefficient. For 
amplitudes A larger than about 3.0, the diffusion coefficients are larger than their 
bare diffusion coefficient values[7,8]. Thus this peak could be attributed to optimised 
enhancement of the escape rate by modulation for a given noise strength[9,10,ll]. It will 
be interesting to see how the peaking behaviour of D gets affected by the variation of 
temperature and frequency. For high frequency drive (ou = 10, for example), however, 
the peak disappears completely (Inset, fig. 5). Not only the peak disappears, the overall 
diffusivity itself gets reduced by about two orders of magnitude indicating that the 
particle excursion gets confined to a single well with only occasional jump to the right 
and left of the barriers, exhibiting same behaviour for both ±V(x). In fig.6, we plot the 
diffusion coefficient for the asymmetric potential V(x) = cos a; + cos 2x + 0.25 sin 3x as 
a function of amplitude A of the ac forcing A cos(o>t) with uj — 1.0 at temperature T = 



0.25. In this case too D shows similar behavior to the symmetric potential case, fig.5. 
Remarkably, also for the asymmetric potential D + (A) = D_(A). 



IV. CONCLUSION 

In conclusion, we have shown that the current and diffusion coefficients, in static 
potentials ±V(x) in the presence of constant force, are related by symmetry relations 
J + (F) = -J_(-F) and D+{F) = D-(-F), respectively. When the Fokker-Planck 
equation of motion is mapped on to the Schrodinger equation these potentials lead 
to supersymmetric Schrodinger partner potentials, whose eigenspectra are isospec- 
tral. Also, eigenfunctions are in one to one correspondence. Consequently, in the 
large time asymptotic limit these differing potentials give the stated behaviour of the 
transport coefficients. We have also shown, numerically, that in the presence of ac 
drive the diffusion coefficients exhibit same behaviour for ±V(x). This is surprising 
because supersymmetry arguments do not apply in the presence of time dependent 
fields. This , therefore, requires further investigation. 
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FIG. 1: The potentials V(x) = cos x + cos 2x (solid line) and —V(x) (dashed line) are plotted 
in dimensionless units for comparision. 




FIG. 2: The potentials V(x) = cos x + cos 2x + 0.25 sin 3a; (solid line) and —V(x) (dashed line) 
are plotted in dimensionless units for comparision. 




FIG. 3: The currents J± and (inset) D± are plotted as a function of constant applied force F 
for symmetric potentials. The subscripts =L refer to the potentials ±V(cc) and at temperature 
T = 0.1. Results for ±V(x) are identical. 




FIG. 4: The currents J± and (inset) D± are plotted as a function of constant applied force 
F for asymmetric potentials. The subscripts ± refers to the potentials ±V(x) and at tem- 
perature T = 0.1. Continuous lines correspond to +V(x). 




FIG. 5: The diffusion constants D± are plotted as a function of the amplitude A of ac force 
for the symmetric potential ±V(x) at frequencies u = 1.0 and (inset) lo = 10.0 and at 
temperature T = 0.25. The continuous line corresponds to D + . 




FIG. 6: The variation of D± as a function of amplitude A of a ac drive for the asymmetric po- 
tential ±V(x), (eq.(14)) at temperature T = 0.25 and high frequency uj = 1.0. Continuous 
line for +V{x). 



